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I. INTRODUCTION 



. Recently quantum teleportation, quan 



QQl 



have spurred 



In the last decade quantum entangled states have showed remarkable applications and 
become one of the key resources in the rapidly expanding field of quantum information 
processing. The history can be traced back to the earlier well-known papers of Einstein, 
Podolsky and Rosen 1], Schrodinger P| and Bell ^. 
tum cryptography, quantum dense coding and parallel computation 
a fiurry of activity in the effort to fully exploit the potential of quantum entanglement. 
Despite of its importance, we do not yet have a full understanding of the physical character 
and mathematical structure for entangled states. We even do not know completely wether a 
generic quantum state is entangled, and how much entanglement remained after some noisy 
quantum processes. 

A state of a composite quantum system is said to be disentangled or separable if it can 
be prepared in a local" or classical" way. A separable bipartite system can be prepared 
as an ensemble realization of pure product states iV'i)^ = 1, M for some positive 

integer M) occurring with a certain probability pf. 

i 

where pf = pf = \4>i)B {4>i\, ^ \^i)B normalized pure states 

of the subsystems A and B, respectively iq. If no convex linear combination exists for a 
given pab, the state is called entangled" and includes quantum correlation. 

For a pure state pab, it is straightforward to judge its separability: a pure state pab is 
separable if and only if there is only one item in Eq. (^Q) and p"^ resp. p^ are the reduced 
density matrices defined as p^ = Ttb^Pab) and p^ = Tta^Pab)- However, for a generic 
mixed state pab-, finding a decomposition like in Eg. (ITl) or proving that it does not exist is a 
non-trivial task (we refer to recent good reviews Isllgllioll and references therein). There has 
been considerable effort in recent years to analyze the separability and quantitative character 
of quantum entanglement. The Bell inequalities satisfied by a separable system give the first 
necessary condition for separability [sl . Many years after the appearance of Bell inequalities, 
Peres made an important step forward in 1996 by showing that partial transpositions with 
respect to one and more subsystems of the density matrix for a separable state are positive. 
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p^- > 0, (2) 

where X is either A or B, p^^ stands for the partial transpose with respect to X. Thus p^^ 
should have non-negative eigenvalues (this is known as the PPT criterion or Peres-Horodecki 
criterion) which was further shown by Horodecki et al. to be sufficient for 2 x 2 
and 2x3 bipartite systems. Meanwhile, these authors also found a necessary and sufficient 
condition for separability by establishing a close connection between positive map theory and 
separability In view of the quantitative character for entanglement, Wootters succeeded 
in computing the ^''entanglement of formation'^ 13] and thus obtained a separability criterion 
'or 2 X 2 mixtures The ^^reduction criterion" proposed independently in ^| and 

gives another necessary criterion which is equivalent to the PPT criterion for 2 x n 
composite systems but is generally weaker. Pittenger et al. gave also a sufficient criterion 
for separability connected with the Fourier representations of density matrices [l^. Later, 
Nielsen et al. j3| presented another necessary criterion called the majorization criterion: 
the decreasingly ordered vector of the eigenvalues for pab is majorized by that of p a or pb 
alone for a separable state. A new method of detecting entanglement called entanglement 
^nesses was give,. i„ H a„d QQ. So,„e necessary a„d sufflcie,. cdter.a of separabUity 

the density matrix are also known [21|, |23, 12^. In addition, it was 



for low rank cases o: 

shown in [2^ and [25I that a necessary and sufficient separability criterion is also equivalent 
to certain sets of equations. A PPT extension based on semidefinite programs is proposed 
in 2^ which can test numerically the separability. 

However, despite these advances, practical and easily computable criteria for a generic 
bipartite system are mainly limited to several ones such as the PPT, reduction, majorization 
criteria as well as the PPT extension. Very recently Rudolph |22| and K. Chen and L.A. Wu 
128^ proposed a new operational criterion for separability: the realignment criterion (named 



given in Ref. 



following the s ugg estion of 29[, it coincides with the computational cross norm criterion 



27|). The criterion is very simple to apply and shows dramatic ability to 



detect many of the bound entangled states j2^ and even genuinely tripartite entanglement 
j^l • Soon after the appearance of , Horodecki et al. showed that the PPT criterion and 



realignment criterion are equivalent to the permutations of the density matrix's indices 
A simple .ngle tpnewoA for the=e cnteHa for the .nuMpartite ca.e m any dimensions was 

recently given in j3fl(| and is called the generalized partial transposition criterion {GPT) which 



includes, as special cases, the Peres-Horodecki criterion (PPT), the realignment criterion 
and the permutation indices criterion for density matrix. Some further properties of the 



realignment criterion have been very recently derived in j31|. 

In this paper we present a systematic generalization of the reduction criterion employing 
a realignment technique of a certain matrix constructed from the density matrix. This 
criterion includes the reduction criterion and the GPT criterion as special cases. It unifies 
them in a single simple framework. Thus our criterion is a strong separability test for 
a generic bipartite or even for multipartite quantum states in arbitrary dimensions. The 
detailed constructions are given in Section ITTl where the reduction and the GPT criteria are 
shown to be two special cases of our new criterion. Some interesting examples are given in 
Section lim A brief summary and some discussions are given in the last section. 

II. THE CRITERIA FOR SEPARABILITY 

In this section we study the separability of the density matrix for any bipartite system in 
arbitrary finite dimension. Motivated by the reduction criterion and the Kronecker product 



approximation technique j32,|33[, we give a new separability criterion by analyzing the trace 
norm for some realigned version of a matrix constructed from the whole density matrix and 
its reduced ones. 

A. Some notation 



We first introduce some notations used in various matrix operations (see e.g., 34l l35l|): 
Definition: For any mx n matrix A = [aij], with entries aij, we define a vector Vvec{A) by 

where t represents the standard transposition operation. Let %. (resp. 7^) denote the row 
transposition (resp. column transposition) of A: 

%.{A) = (V.ec(A))*, (3) 

Z{A)=yvec{A). (4) 



It is easily verified that 



%%.{A) = %%{A) = (5) 



For example, for a 2 x 2 matrix A, 



A 



ail ^12 



0-21 (^22 



we have 



Oil 021 



an 



^21 



I3'12 I3'22 



(6) 



(7) 



ai2 

y <3'22 y 

For a generic matrix A = ^^Aj K) (j| = ^^Aj \i) ® (j| = ^^Aj (j| ® where \i) , \ j) 

i,j id i,j 

are vectors of a suitably selected normalized orthogonal basis, {i \ , (j| are the corresponding 
transpositions (not conjugate transpositions). Under the operations %. and %. one has: 



A 5^ A,, (j| ® (z| ^ |j) ® (z| = A*, 

A ^ 5^ A,, |j) ® ^ 5^ A,, |j) ® (^1 = A*. 



(8) 

_ (9) 

We further define (resp. Xi^) {k = A, B) to be the row (resp. column) transposition with 
respect to the subsystems A,B. We set T^xi,x2,...} = %:i%c2--- ^i; ^2 £ {'^a, c^i, r^, c^}. A 
generalized partial transposition operation [GPT operation) for a bipartite density matrix 
is thus given by 7y, y C {rA^CA^rB^CB}, where Ty stands for all partial transpositions 
contained in y which is a subset of {r^, c^i, r^, c^}. With these notations, the realignment 
criterion can easily be stated. For example, for a bipartite system, we only need to make 
partial transpositions with respect to 3^ = {ca.tb}- This is equivalent to the reahgnment 
operation given in j2^ and j28|, for the proof, see jS^. 

We also need the following result in matrix analysis. Let Z be an m x m block matrix 
with block size nxn. We define a realigned matrix Z of size x that contains the same 
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elements as in Z but in different positions, 



Vt,ec(^l,l)* 
^vec(-^m,l) 
"^vec ( ^1 ,m ) 

Z has the singular value decompositions: 

z = uj:v^ = ^2"" 



(10) 



(11) 



i=l 



where U = [uiU2 ■ ■ -Um^] € C"^^^™"' and V = [viV2 ■ ■ -fna] G (C"^^"^ are unitary matrices, S 
is a diagonal matrix with elements cri > (^2 > ■ ■ ■ > Cg > and q = min(m^,n^). In fact, 
the number of nonzero singular values cXj is the rank r of the matrix Z, and are exactly 
the nonnegative square roots of the eigenvalues of ZZ'' or Z'^Z [35]. Based on the above 
constructions, Z can be expressed as: 



z = Y,ix^^yi} 



(12) 



with Vvec{Xi) = ^/a'iUi and Vvec{Yi) = ^falvl 



i=l 

HQ- 



B. The generalized reduction criterion for separability 

We now derive the main result of this paper: a generalized reduction criterion for sepa- 
rability of bipartite quantum systems in arbitrary dimensions. 



1. The main theorem 



The reduction criterion given in 







15l | and |l6| says that for any bipartite m x n separable 



states, the following inequalities should be satisfied simultaneously: 



Im ® - Pab ^0, Pa 'S)ln- Pab > 0, 



(13) 
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where pa,b are the reduced density matrixes with respect to the subsystems A and S, Im 
(resp. In) is an m (resp. n) dimensional identity matrix. This criterion is shown to be 
equivalent to the PPT criterion for 2 x n system but it is generally weaker than the PPT 



criterion 
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So it certainly cannot detect any bound entangled states which are PPT. 
Noticing this fact and the powerful ability of the realignment criterion, in particular its 
generalization: the GPT criterion, we expect that some stronger test may appear. The 



reduction criterion is in essence a positive map. Combining the technique 



map and the Kronecker product approximation technique for a matrix 
more general map: 

PAB ^ PAB = abimn " Cl^m ® Pb - bpA ® !„ + PAB, 



or this positive 
3^, we apply a 

(14) 



where a, b are arbitrary complex numbers. We are going to derive a necessary separability 
condition of Pab in terms of the trace norm (Ky Fan norm) of a matrix obtained from a GPT 
map on Pab- The trace norm of a matrix Z is a unitary invariant norm which is defined as 
the sum of all the singular values of Z, or alternatively the sum of the nonnegative square 
roots of the eigenvalues of ZZ'' or Z'^Z. We thus arrive at the following separability criterion 
for a bipartite system: 

Theorem: // a bipartite density matrix pab defined on an m x n space is separable, then 
the generalized reduction version pab of pab should satisfy 



\Pab 



< hnh 



Wy C {rA,CA,rB,CB}, 



(15) 



where 7^^, or 7^^ (k = A, B) stands for transpositions with respect to the row or column for 
the subsystem k. The numbers ha,hf, are defined by 



|a - 1| + (m - l)|a|, r^, G 3^ or rA,CA^y, 
(|a - Ip + (m - l)|ap)5, G 3^, Ca ^ y or 

CA ey, TA^ y, 



hh = 



1 



\b - 1\ + {n - l)\b\, rB,CBeyor rB,CB^y, 
i\b - Ip + {n- l)|6p)5, rBey, CB^y or 

cb ey, tb^ y, 



where a, b G C, Ty represents partial transpositions with respect to every element contained 
in the set y of the related subsystems. 

Proof: We only need to find the bound for the trace norm of \\pab ^\\ with respect to 
some operations Ty for any separable states. Considering a separable bipartite system, we 
suppose pab has a decomposition, pab = J2iPiPf ® pf with < < 1, XltP* = 1- Under 
map (jl4j) it is evident that 

Pab Pab = ablmn - a^m ® Pb - bp a ®ln + Pab 

= J2iPMm- Pt)®ibln- pf), 

where pA and ps are the reduced density matrixes defined by pA = TrsipAB) = T^iPipti 
pB = TvAipAB) = Y^iPiPi- We have 



Pi ® P. 
Ui 



where we have diagonalized the (rank one) density matrix pf (resp. pf) with the m (resp. 
?T,)-dimensional unitary matrix Ui (resp. Vj). It is straightforward to check that 

(al™ - pf) ® {bin - pf) = U.AUl ® V.BVf 







V 



'i 1 



(16) 



where A and B are diagonal matrices: 

( . 1 



A 



B 



For an m X m (resp. n x n) matrix P (resp. Q) acting on the complex space associated 
to the subsystem A (resp. B), the trace norm of the matrices P and Q has the property: 
||-P®<5|| = ||-P|| ■ IIQII- And for any matrices X, Y Z acting on the subsystem A (or B), 
the Vyec operations have the property: V^ed^y Z) = {Z'^ ® X)Vvec{y) , where both sides 
are column vectors and the tensor operation ® has nothing to do with different subspaces 
A and B b 
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Let pab ^ denote the transformed matrix of pab under the partial transposition Ty. 
Without loss of generality we suppose that we only make a row transposition to the subsys- 
tem A. According to Q, we have 

= {V,eMAU})Y®ViBV^ 

= {Vvec{A))\U,^ ® Ul) ® V^BVl (17) 
In terms of the unitary invariant property of the trace norm we obtain 

\mAU} (^ViBV^fi^A}\\ 

= ||(V.ec(^))*(f/.t®f/*)||.||W/|| 

= \\{y,,M)n-m\ 

= (|a - + (m - l)|a|2)^(|6 -l\ + {n- l)\h\). 
Using the convex property of the trace norm, we get 

= 1 1 Y^Pi^UiAU} ® ViBV^fi^A} 1 1 

i 

< iJ2P^)i\a - + (m - lM)H\b -l\ + in- l)\b\) 

i 

= (|a - + (m - l)|ar)^(|6 - 1| + (n - 1)|6|). (18) 

A corresponding procedure can be applied for the column transposition to the subsystem A 
and corresponding operations for the subsystem B: 

< (|a - 1|2 + (m - l)|a|2)l • (|6 - 1| + (n - 1)|6|), 

< {\h - 1|2 + {n- l)\h\^)^ ■ (|a - 1| + (m - l)|a|), 

< (|6 - 1|2 + (n - l)|6nl . (|a - l| + (m - l)|a|). 

For the operations with respect to 3^ = {r^i, c^}, we have in fact the PPT operation to the 



subsystem A of pab and thus 



= \\J2PiiU.AUjy^V,BV^\\ 

i 

i 



= (|a-l| + (m-l)|a|)(|6-l| + (n-l)|6|). (19) 
If the subsystem A is left unchanged, i.e., both ca,ta ^ y, we have 

IIpabII 

= \ \Y,P^U^AUj0ViBV^\\ 

i 

< ^p^\\A® B\\ 

i 

= (|a-l| + (m-l)|a|)(|6-l| + (n-l)|6|). (20) 

For any other combinations of y from {r^, ca, tb, cb}, following the same procedure above, 
we arrive at the final result p5|) . ■ 

2. Special cases reducing to other necessary criteria 

We show now that the Theorem actually encompasses two previous strong computational 
criteria for separability. 

a. The reduction criterion For the case of a = and 6 = 1, we have 

M— Tvi, ^ I ('^-l)' rB,CBey OTrB,CB^y, , , 

Wpab "^11 < < , (21) 

1^ (n - 1) 2 , rB e y or Cb e y. 

When rB, cb ^ 3^, we have further 

Wpa^K - PabW < n - 1. (22) 
For the case of a = 1 and 6 = 0, one obtains similarly 

W^m^ Pb - PabW <m - 1. (23) 
10 



Eqs. fl22jl and (j23|l are equivalent to the reduction criterion, since for separable states pos- 
itivity of Eq. p3|) means that the trace norm is the sum of the eigenvalues, that is, the 



singular values, due to the Hermitian property of pa ® \ ~fiAB and Im ® — Pab [35 |. 

h. The GPT criterion The GPT criterion derived in jSOf says that for any GPT op- 
erations, the trace norm for the realigned matrix is not greater than 1. Violation of that 
inequality means existence of entanglement. The GPT criterion includes the realignment 
criterion and the PPT criterion as special cases. Now the GPT criterion is just one of the 
special cases of our generalized reduced criterion: a = and 6 = 0. In this case we have 

||p5^||<1, V3^C W,CA,rB,CB}. (24) 
This is exactly the GPT criterion. 



3. Invariance of our generalized reduction criterion under local unitary transformations 

The trace norm ||pab'^"'^|| is invariant under local unitary transformations. To see this, 
let Wa_ (resp. Wb) be unitary transformations on the subsystem A (resp. B). Under the 
local unitary transformation Wa®Wb, Pab is mapped to = {Wa^Wb)pab(}V\®W^). 
If for any GPT operations, the local unitary transformation only contributes some unitary 
factors to pab ^ , "we would certainly have | \pab 1 1 = I \Pab II; unitary invariance 

of the trace norm. In fact, pab has the decomposition pab = XlLi ® where q = 
min(m^,n^), in terms of the procedure (fTUI) to (fT^ . Without loss of generality, we consider 
a row transposition on the subsystem A. Setting 7j = {Wa ® WB){o:i (8> l3i)(W\ ® W^), we 
have 

7?"^' = V.ec{WAa,W\y ® WBml 
= {V.ec{a^)y{WA^ ® W'a) ® WBml, 

= WB{{V,ec{ai)y ^ f3,){WA^ ^W'^^Wl) (25) 

Summing over all the components 7j '^'^ , i = 1,2, we have p\^ = 

WbPab'^^'^^^ {Wa^ ® W\ ® Wb)- Therefore p\^ ^ and Pab'^^''^^ are equivalent under the 
unitary factors Wb and {Wa'' (S> W^ (S> Wj^), which keep the trace norm invariant. 

The same procedure can be used to perform column transposition, partial transposition 
of some subsystems, and any combinations of these GPT operations, to show that the trace 
norm is an invariant under local unitary transformation. 
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III. EXAMPLES 



In this section, we consider two examples to illustrate some special characters of the 
criterion compared with the previously known reduction criterion and the GPT criterion. 

Example 1: 3x3 Werner state. 

Consider the family of ci- dimensional Werner states [3], 

where — 1 < / < 1, V{a ® /3) = /3 ^ a, the operator V has the representation V = 
Ylt'jLo cind Wd is non-separable if and only if — 1 < / < 0. As is well known, the 

entanglement in a 2 x 2 Werner state can be detected completely using PPT, reduction 
and realignment criteria. But for higher dimensions the reduction criterion fails while PPT 

2 



succeeds. The realignment criterion can recognize entanglement for — 1 < f < 4 — 1 j27|. 
Here for simplicity we consider the 3x3 Werner state given by (j26|) and take a,b & M. We 
plot N = max{||p^B {^a.'-s}!! _ hahh,0} as a function of b and / for the cases of a = and 
a = 1 respectively, FigHJ 

For the case a = 0, we see that > for all — 1 < / < — | and 6 = or |. From the 
Theorem, we have /la/if, = 1 for 6 = or | when a = 0. pab is in fact the same for 6 = 
or I and we obtain = max{ ^^^'^.^"'^ , 0} by direct computation. For that case a = 1, we 
still have > for all — 1 < / < — | but with 6 = — | or 1. Also we see that hahf, = 2 for 
6 = — I or 1 when a = 1. pab is the same for both 6 = — | and 1. In this case we have again 

= max{i^— 0} by direct verification. 

Example 2: Horodecki 3x3 bound entangled state 

n 

Horodecki gives a very interesting weakly entangled state in [36|| which cannot be detected 
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FIG. 1: Depiction of = max{| IpIi'^'^A-'s) 1 1 - hahh,0} for a 3 X 3 Werner state as a function of 
b and / when a = (the top figure) and a = 1 (the bottom figure), respectively. 

by the PPT criterion. The density matrix p is real and symmetric: 
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FIG. 2: Depiction of = max{||pyiB ^'^^-'s} || — /i^j/ij,, 0} for a Horodecki 3x3 bound entangled state 
as a function of h and c when a = (the top figure) and a = 1 (the bottom figure), respectively. 

where < c < 1. The entanglement in this state is very difficult to detect with previous 
operational criterion in general. In [s^, Horodecki showed that the range criterion could 
recognize the entanglement. For this state, the simple realignment criterion and its general- 
ization: the GPT criterion can, surprisingly, detect the entanglement for all permissable c, as 
shown in 2^. Here, we will observe the behavior in the language of the generalized reduction 
criterion (for real values of a and 6, as in Example 1). The partial transposition pab '^''^'"^^ 
using our Theorem fails to detect the entanglement by direct calculation. Therefore we are 
only concerned with the case of pab '''-^'''s^ which corresponds to the realignment of Pab- 
For the case of a = and a = 1, respectively, we plot = max{||p^B {'=a.'-b}|| — hahb,0} as 
a function of b and c in Fig|21 For the case of a = 0, we see that > for all < c < 1 and 
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6 = or |. In this case hahb = 1 for 6 = or | when a — 0. But pab has different forms for 
6 = and |. For the case a — 1 we still have A^" > for all < c < 1 but with 6 = — | or 1. 
Also we see that hahb — 2 and pab has quite different forms for 6 = — | and 1 when a = 1. 

However, the function N has the same value for the case a = when 6 = or |, or the case 
a = 1 when b = — | or 1. 

The results of the above two examples arc a little bit surprising compared with the 
relationship between the PPT and reduction criteria. As is well known, the PPT criterion 
is equivalent to the reduction criterion for a 2 x n system. That is to say, for a = and 
6 = 1, or a = 1 and 6 = they give the same result for a 2 x n system. But in the case of 
the 3x3 system, we have identical results for different b other than the value of 0, 1 if we 
apply the realignment operations to Pab- This interesting phenomenon also occurs in higher 
dimensional cases. So our generalized reduction includes all the results from the original 
realignment or GPT criteria but recognizes entanglement in a very subtle way which is quite 
different from other criteria. 

IV. CONCLUSION 

Summarizing, we have introduced a computational criterion, which we call the "gen- 
eralized reduction criterion", providing a necessary condition for separability of bipartite 
quantum systems in arbitrary dimensions . This criterion combines many virtues of the 
reduction criterion and the GPT criterion. It gives a unified framework for the two criteria 
and provides a powerful necessary condition for separability using just simple matrix opera- 
tions. Some interesting characters of this criterion are showed by two typical examples. The 
theorem can be straightforwardly generalized to the multipartite case by introducing more 
free parameters like a and b in the theorem. 

We expect that this construction not only expands theoretically our sight in detecting 
entanglement of a general quantum state, but also sheds some light on possible ways to the 
final solution of the separability problem. We conjecture that any future stronger operational 
separability test should in principle include as special cases the simple GPT criterion, in 
particular the PPT criterion which is necessary and sufficient one for 2 x 2 and 2x3 system. 
This paper is an attempt following this way, though we have not yet found an example of a 



15 



bound entangled state which can be detected by this criterion but not by the GPT criterion. 
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